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Abstract. We study the behavior of max-algebraic powers of a reducible 
nonnegative matrix A S R^ xn . We show that for t > 3n 2 , the powers A t 
can be expanded in max-algebraic sums of terms of the form CS t R, where C 
and R are extracted from columns and rows of certain Kleene stars, and S is 
diagonally similar to a Boolean matrix. We study the properties of individual 
terms and show that all terms, for a given t > 3n 2 , can be found in 0(n 4 logn.) 
operations. We show that the powers have a well-defined ultimate behavior, 
where certain terms are totally or partially suppressed, thus leading to ulti- 
mate CS t R terms and the corresponding ultimate expansion. We apply this 
expansion to the question whether {A l y, t > 0} is ultimately linear periodic 
for each starting vector y, showing that this question can be also answered in 
0(n 4 logn) time. We give examples illustrating our main results. 



1. Introduction 

By max algebra we understand the analogue of linear algebra developed over the 
max-times semiring K max ,x which is the set of nonnegative numbers R+ equipped 
with the operations of "addition" a<S)b :— max(a, b) and the ordinary multiplication 
a®6:=ax!). Zero and unity of this semiring coincide with the usual and 1. The 
operations of the semiring are extended to the nonnegative matrices and vectors in 
the same way as in conventional linear algebra. That is if A — (ay), B = (bij) and 
C = (cij) are matrices of compatible sizes with entries from R + , we write C = A(BB 
if Cij = dij © bij for all i,j and C = A <g> B if = fc a lk b kj = max. k (a ik b k j) for 
all i,j. 

If A is a square matrix over R + then the iterated product A<E> A®...® A in which 
the symbol A appears k times will be denoted by A . These are the max-algebraic 
powers of nonnegative matrices, the main object of our study. 

The max-plus semiring M m ax,+ = (1U {— oo}, = max, ® = +), developed over 
the set of real numbers R with adjoined element — oo and the ordinary addition 
playing the role of multiplication, is another isomorphic "realization" of max alge- 
bra. In particular, x exp(a;) yields an isomorphism between R maXj+ and R m ax,x- 
In the max-plus setting, the zero element is — oo and the unity is 0. 

The main results of this paper are formulated in the max-times setting, since 
some important connections with nonnegative and Boolean matrices are more trans- 
parent there. However, the max-plus setting is left for the examples in the last 
section, in order to appeal to the readers who work with max-plus algebra and 
applications in scheduling problems and discrete event systems [TJ [17] . 

The Cyclicity Theorem is a classical result of max algebra, in the max-times set- 
ting it means that the max-algebraic powers of any irreducible nonnegative matrix 
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are ultimately periodic (up to a scalar multiple) , with the period equal to cyclicity 
of the critical graph. This theorem can be found in Heidergott et al. [I7| Theorem 
3.9], see also Baccelli et al. [TJ Theorem 3.109] and Cuninghame-Green [5J Theorem 
27-6] (all stated in the max-plus setting). However, the length of the pre-periodic 
part can be arbitrarily large and the result does not have an evident extension to 
reducible matrices. 

The behavior of max-algebraic powers of reducible matrices relies on connec- 
tions between their strongly connected components and the hierarchy of their max- 
algebraic eigenvalues. This has been well described in a monograph of Gavalec |16j . 
see also QUID]- 

The relation of the Cyclicity Theorem to the periodicity of Boolean matrices is 
understood but not commonly and explicitly used in max algebra. For instance the 
construction of cyclic classes [21 E] appears in the proof of Lemma 3.3 in [17] (without 
references to literature on Boolean matrices). This relation becomes particularly 
apparent after application of a certain D~ 1 AD similarity scaling called visualization 
in [25j [26] , see also [TO] Qj] . Semancikova [22j [23] realized that the cyclic classes of 
Boolean matrices are helpful in treating computational complexity of the periodicity 
problems in max-min algebra, with analogous max-algebraic applications in mind 

m- 

A result by Nachtigall ;2Tj states that, though the length of the preperiodic part 
cannot be polynomially bounded, the behavior of matrix powers after 0(n 2 ) can 
be represented as max-algebraic sums of matrices from certain periodic sequences. 
Molnarova [19] studies this Nachtigall expansion further, showing that for a given 
matrix power after 0(n 2 ) the representing Nachtigall matrices can be computed in 
0(n 5 ) time. 

In the general reducible case, the sequences of entries {a\j,t > 0} of A* are 
ultimately generalized periodic 9. 20 . meaning that for t > T where T is sufficiently 
large, they may consist of several ultimately periodic subsequences which grow 
with different rates. Gavalec [T5J [TB] showed that deciding whether {a\^,t > 0} is 
ultimately periodic is in general NP-hard. In a related study, Butkovic et al. [5] 
[6] considered robust matrices, such that for any given x the sequence is 
ultimately periodic with period 1. The conditions formulated in [6] can be verified 
in polynomial time, which suggests that such "global" periodicity questions must 
be tractable. 

The main goal of this paper is to find a common ground for the above pieces 
of knowledge on matrix periodicity in max algebra. We introduce the concept of 
CSR expansion, in which a matrix power A' is represented as the max-algebraic 
sum of terms of the form CS t R, called CSR products. Here C and R have been 
extracted from columns and rows of certain Kleene star (the same for both), and 
S can be made Boolean by a similarity scaling D~ 1 SD. The matrix CR appeared 
previously as spectral projector [U[7], and S typically arises as the incidence matrix 
of a critical graph. 

After giving necessary preliminaries in Section (2] we start in Section [3] by study- 
ing the CSR products. We show that they form a cyclic group and describe the 
action of this group on the underlying Kleene star. We also emphasize the path 
sense of these operators, see Theorem 13.31 thus providing connection to the ap- 
proach of QH OH [HI [22] [23]. In Section H we establish the algebraic form of 
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the Nachtigall expansion which controls the powers A t after t ~ 0(n 2 ). See The- 
orem 14.21 In Section [5] we show that at large t certain Nachtigall terms become 
totally or partially suppressed by heavier ones, leading to the ultimate expansion of 
matrix powers. This result, see Theorem l5.6[ can be understood as a generalization 
of the Cyclicity Theorem to reducible case. In Section[6]we treat the computational 
complexity of computing the terms of CSR expansions for a given matrix power, 
showing that in general this can be done in no more than 0(n 4 log n) operations. In 
Section[7J which extends the results of Butkovic et al. 5, 6 on robust matrices, we 
describe orbit periodic matrices A, i.e., such that the orbit A y is ultimately linear 
periodic for all initial vectors y, see Theorem 17.61 We use the ultimate expansion 
to show that the conditions for orbit periodicity can be verified in no more than 
0(n 4 log n) operations, see Theorem 17.81 and its Corollary. We conclude by Section 
E] which contains some examples given in the max-plus setting. 

2. Preliminaries 

In this section we recall some important notions of max algebra. These are the 
maximum cycle geometric mean, the critical graph and the Kleene star. We close 
the section with nonnegative similarity scalings and cyclic classes of the critical 
graph. 

Let A = (aij) e M" x ". The weighted digraph V{A) = (N(A),E{A)), with the 
set of nodes N(A) — {l,...,n} and the set of edges E(A) — | ay ^ 0} 

with weights w(i,j) — ay, is called the digraph associated with A. Suppose that 
P = {i\,...,i p ) is a path in Da, then the weight of P is defined to be w(P) — 
■ • • &i p _ 1 i p if p > 1, and 1 if p = 1. If i\ = i p then P is called a cycle. The 
length of P, denoted by l(P), is the number of edges in P (it equals p—1 here). 

When any two nodes in D(A) can be connected to each other by paths, the matrix 
A is called irreducible. Otherwise, it is called reducible. In the reducible case, there 
are some (maximal) strongly connected components of D(A), and a number of 
nodes that do not belong to any cycle. We will refer to such nodes as to trivial 
components of T>(A). 

The maximum cycle geometric mean of A, further denoted by X(A), is defined 
by the formula 

(2.1) X(A) =max( W (F c )) 1 / fe , 

where the maximization is taken over all cycles P c = . . . , ik), for fc = 1, . . . , n, 
in the digraph T>{A). 

The Cyclicity Theorem (JTT1 Theorem 3.9], see also Q] 0) states that if A 
is irreducible then after a certain time T(A), there exists a number 7 such that 
A t+1 = y(A)A t for all t > T(A). Thus A is the ultimate growth rate of matrix 
powers in this case. If A = 1 then A t+1 = A* for t > T(A), in which case we say 
that {A*, t > 0} is ultimately periodic. 

Remarkably X(A) is also the largest max-algebraic eigenvalue of A, meaning the 
largest number A for which there exists a nonzero x S M™ such that A (g> x = Xx, 
see [2 IH [17] and references therein. 

The operation of taking the maximal cycle geometric mean (m.c.g.m. for short) 
is homogeneous: X(aA) — aX(A). Hence any matrix, which has X(A) ^ meaning 
that V(A) is not acyclic, can be scaled so that X(A/X(A)) = 1. Following [3], 
matrix A g IR™ xn with X(A) = 1 will be called definite. 
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A cycle P c = . . . , i k ) in T>(A) is called critical, if (w(P c )) 1/k = X(A). Every 
node and edge that belongs to a critical cycle is called critical. The set of critical 
nodes is denoted by N C (A), the set of critical edges is denoted by E C (A). The 
critical digraph of A, further denoted by C(A) = (N C (A), E C (A)), is the digraph 
which consists of all critical nodes and critical edges of T>(A). 

The cyclicity of an irreducible graph is defined as the g.c.d. (greatest common 
divisor) of the lengths of all its simple cycles. The critical graph defined above may 
have several strongly connected components, and in this case the cyclicity is the 
l.c.m. (least common multiple) of their cyclicities. This gives the number 7 which 
appears in the Cyclicity Theorem [TJ [17] , and it can be shown that the ultimate 
period cannot be less than 7 (in particular, this follows from the approach of the 
present paper). 

There is no obvious subtraction in max algebra, however we have an analogue 
of (I - A)- 1 defined by 



where / is the identity matrix. This series converges to a finite matrix if and only 
if \{A) < 1 [TJ IH [17], and then this matrix A* = (a*-) is called the Kleene star 
of A. This matrix has properties (A*) 2 = A* and, clearly, A* > I. It is important 
that the entries of max-algebraic powers A k — (a^ ) express the maximal weights 
of certain paths: a k j is equal to the greatest weight of paths P that connect i to j 
and have length k. The entry a*j for i 5^ j is equal to the greatest weight of paths 
that connect i to j with no restriction on their lengths. 

As in the nonnegative linear algebra, we have only few invertible matrices, in 
the sense of the existence of (nonnegative) A^ 1 such that A^ 1 <E> A — A <E> A^ 1 = I. 
More precisely, such matrices can be diagonal matrices 



for a positive x = (x±, . . . ,x n ), or monomial matrices obtained from the diagonal 
matrices by permutations of rows or columns. Nevertheless, such matrices give 
rise to very convenient diagonal similarity scalings A 1— > X~ l AX. Such transfor- 
mations do not change A (A) and C(A) [T2]. They commute with max-algebraic 
multiplication of matrices and hence with the operation of taking the Kleene star. 
Geometrically, they correspond to automorphisms of K™ , both in the case of max 
algebra and in the case of nonnegative linear algebra. The importance of such 
scalings in max algebra was emphasized already in J3 , Ch. 28. 

By an observation of Fiedler and Ptak 13J, for any definite matrix A there is 
a scaling X such that X~ x AX is visualized meaning that all critical entries of the 
matrix equal 1 and all the rest are less than or equal to 1. It is also possible to 
make X~ l AX strictly visualised [26 , meaning that only critical entries are equal 
to 1. If a matrix is visualised, or strictly visualised, the same is true for all powers 
of this matrix, meaning that the critical graph can be seen as a Boolean matrix 
that "lives by itself" . Thus there is a clear connection to the powers of Boolean 
matrices. 



(2.2) 



A* := I® A® A 2 ® 



(2.3) 
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The periodicity of powers of Boolean matrices is ruled by cyclic classes [2 also 
known as imprimitivity sets [3], which we explain below. We note that this notion 
appeared already in a work of Frobenius [14] . 

Proposition 2.1 (e.g. Brualdi-Ryser [3 j). Let G — (N,E) be a strongly connected 
digraph with cyclicity jg ■ Then the lengths of any two paths connecting i G N to 
j G N (with i,j fixed) are congruent modulo jg- 

Proposition ^, ll implies that the following equivalence relation can be defined: i ~ 
j if there exists a path P from i to j such that l(P) = (mod 7c). The equivalence 
classes of G with respect to this relation are called cyclic classes [H [22l [23] . The 
cyclic class of i will be denoted by [i] . 

Consider the following access relations between cyclic classes: [i] —>t [j] if there 
exists a path P from a node in [i] to a node in [j] such that l(P) = t(mod7c). In 
this path P with l(P) = i(mod7c) exists between any node in [i] and any 

node in [j). Further, by Proposition 1 2 . 1 1 the length of any path between a node in 
[i] and a node in [j] is congruent to t, so the relation [i] — » t [j] is well-defined. 

Cyclic classes can be computed in 0(|£?|) time by Balcer-Veinott digraph con- 
densation [2], where \E\ denotes the number of edges in G. At each step of this 
algorithm, we look for all edges which issue from a certain node i, and condense all 
end nodes of these edges into a single node. Another efficient algorithm is described 

Let S = (sij) be the incidence matrix of G, meaning that Sij = 1 if G E 
and Sij — otherwise. 

The ultimate periodicity of such Boolean matrices has been well studied. If 
jG = 1 then the periodicity of 5* starts latest after the Wielandt number W(n) := 
(n - l) 2 + 1 [T8]. This bound is sharp and is due to Wielandt [28] . 

If 7 > 1 then there are even better sharp bounds due to Schwartz [57j. Assume 
w.l.o.g. that S is irreducible, and let n — cry + 1. If a > 1 then the periodicity 
starts at most after W(a)-f + 1 which does not exceed ^- + 7. If a — 1 then it starts 
almost "straightaway", after at most max(l,f). 

3. CSR PRODUCTS 

In this section, given a nonnegative matrix A G R™ xn , we consider max- algebraic 
products of the form CS t R, where S is associated with some subdigraph of the 
critical graph C(A), and matrices C and R are extracted from a certain Kleene star 
related to A. We show that these products form a cyclic group and study their 
periodic properties. We also show that CS t R are related to a distinguished set of 
paths which we call C-heavy. 

We start with a remark that the concept of cyclic classes discussed in Section [2] 
can be generalized to completely reducible digraphs, which consist of (possibly sev- 
eral) strongly connected components, not connected with each other. Importantly, 
the critical digraph of any A G M" x " with A (A) > is completely reducible. 

Let A G R" x ™ have A (A) = 1. Consider any completely reducible subdigraph 
C = (N c , E c ) of C(A). In particular, C may consist of several disjoint cycles of C(A), 
or we can take a component of C(A), or we may just have C = C(A). Denote by 
7 the cyclicity of C and take B := (A 1 )*. Define the matrices C = (cy) G E" xn , 
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FIGURE 1 . The scheme of C and R defined in (|3.1I) 

(sy) G K 1 -" by 



£™ x ™ and S 



0. 



if j G iV c 
otherwise, 

if (ij) G E c 
otherwise. 



0. 



if i e N c 
otherwise, 



The nonzero entries of C, respectively R, can only be in the sub-matrix of B = (A 1 )* 
extracted from columns, respectively rows, in N c . All nonzero entries of S are in the 
principal submatrix Sn c n c extracted from rows and columns in iV c . See Figure Q] 
for a schematic display. 

We show in the next proposition that S can be assumed to be — 1. It can be 
also deduced from the results in [12]. 

Proposition 3.1. Let A G M™ x ™ have X(A) = 1 and let S be defined as above. 
There exists a positive z S K™ such that D~ 1 SD for D :— diag(z) is a — 1 
matrix. 

Proof. As X(S) = 1, we can take z := ©™ =1 5*. This vector is positive, and observe 
that Sz < z since SS* < S*. From this and V(S) = C(S) = C, it can be deduced 
by multiplying z~ 1 SijZj < 1 along cycles that z^ 1 SijZj = 1 for all £ D(S), 



while z i 1 SijZj — for all (i,j) £ T>(S). 



□ 



As S can be scaled to — 1 matrix, we conclude that {S 1 *, t > 0} becomes 
periodic at most after the Wielandt number W(n c ) = (n c — l) 2 + 1, where n c is 
the number of nodes in iV c . Matrix A will be called S -visualized, if S defined in 
(|3.ip is Boolean. In the S'-visualized case, the asymptotic form of 5* for t > T is 
determined by the cyclic classes of C [3]: 



(3.2) 



s t fl, if [*]-*[?'], 
iJ 1 0, otherwise. 

Now we study the Ci'i? products 

(3.3) := C5*i?, 

assuming that A G R+ Xn has A (A) = 1. 
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We start by the observation that if T is the number after which {5*, t > 0} 
becomes periodic, then 

(3.4) S h R = R, CS h = C, V2 7 > T. 

Indeed, (|3.2j) implies that all diagonal entries of S*' 7 with indices in N c are 1 if 
Z 7 > T, which implies S l ~<R > R and CS h > C. On the other hand, S < A and 
so S l ~< < (A~>)*. Further, (A 7 )*P < R and C(A 7 )* < C since (A 7 )*(A 7 )* = (A 7 )*, 
and so S h 'R < (A 7 )*P < R and CS"' 7 < C(A 7 )* < C. 

As {5*, t > T} is periodic, so is {V {t \ t > T}. Moreover, we conclude from 
(|3.4|) that this periodicity starts from the very beginning. 

Proposition 3.2 (Periodicity). p(*+ 7 ) = V {t) for all t > 0. 

Proof. It follows from Eqn. (031) that V {t+h) = V {t) and -pl t +V+ 1 h) = V {t+l) for 
1-1 > T. But p('+ i7 ) = p(t+(i+ih) ) which implies that also 7>( t+7 ) = □ 

It is also useful to understand the meaning of in terms of paths. Given a 
set of paths II, we denote by u>(II) the greatest weight of paths in II, assuming 
w(Jl) = if II is empty. A path will be called C -heavy if it goes through a node in 
N c . The set of all C-heavy paths on T>(A) that connect i to j and have length t will 
be denoted by il' 1 ^ t . We also denote by r the maximal cyclicity of the components 
oiC(A). 

Theorem 3.3 (C-Heavy Paths). Let A G M+ x ™ have X(A) = 1 and let T > be 
such that {S , t > T} is periodic. 

1. Fort> 0, 

(3-5) w(I§ tt ) <V$ . 

2. Fort > T + 2r(n- 1), 

(3.6) w(I$ jtt )>vV . 

Proof. By Proposition 13.11 there exists a diagonal matrix D such that D~ 1 AD is 
S- visualized. As both sides of (I3.5[) and (|3.6[) are stable under similarity scaling of 
A, we will assume that A is already S- visualized. 

1. : Let P e n^ t and w(P) ^ 0. We need to show that 

(3.7) w(P) < V\f 

Path P can be decomposed as P — Pbcg ° -fend, where Pbog connects i to a node 
m G N c and P on d connects m to j. Adjoining to P any sufficiently large number of 
cycles of C that go through m and whose total length is a multiple of 7, we obtain 
a path P' that we can decompose as P' = P^ cg o P( nt o P c ' nd , where P^ eg connects i 
to mi G C and l(Pb cg ) is a multiple of 7, P;' nt connects mi to m 2 G A^, has length 
t and belongs entirely to C, and P c ' nd connects m 2 to j and ^(P e ' nd ) is a multiple of 
7. We conclude that w(Pb cg ) < c imi , w(P c ' nd ) < r m2j - and w(P/ nt ) = s„ im2 = 1. 
We obtain w(P) — w(P') < c imi s t mim2 r m2 j, which implies (|3.7|) . and hence (|3.5j> . 

2. : There exist indices mi and m 2 such that P^ = Ci mi s t mirrl2 r m2 j. This is the 
weight of a path P decomposed as P = Pb og Put Pend, where Pb eg connects i 
to mi, P n t connects mi to m 2 and P on d connects m 2 to j. Here P; n t has length 
t and belongs to the component of C which we denote by T and whose cyclicity 
we denote by %. The lengths of Pb eg and P en d are respectively ^(Pb og ) = and 
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l(Pend) = hir for some h,h, since Z(Pb cg ) and l(P cn d) are multiples of 7 and 7 is 
itself a multiple of n. Paths Pbeg and P en d correspond to certain paths on V(A 7r ) 
with lengths l\ and ^2- If /1 > n or l 2 > n, then we can perform cycle deletion 
(w.r.t. V(A n )) and obtain paths P 2 cg an d -f^nd with lengths kin and k^ where 
ki < n and k 2 < n. For the resulting path P 2 :— P 2 0g ° Pint P c 2 n d we wm 
have w(P 2 ) > w(P) since A = 1. Now we have 2(P b 2 cg ) + l{P^ nd ) < 2r(n - 1). 
If t > T + 2r(n — 1), then the principal submatrix of <!?* corresponding to the 
component T coincides with that of S t ~^ kl+k ^'" , which implies that P„t can be 
replaced by a path Pj 3 t with length t — (fci + fe)"", so that u>(P 3 ) — w(P 2 ) > vf^ 
where P 3 := P b 2 cg o P£ t o P 2 nd G . □ 

This "path sense" of V^' simplifies the proof of the following important law. 

Theorem 3.4 (Group law), p( f i+ f =0 = p(*i)p(fe) f or a llt u t 2 > 

Proof. As (A 7 )*(A 7 )* = {A 1 )*, we have (PC),; = {A 1 )^ = 1 for all i e 7V C , and 
hence PCS' > S. We use this to obtain that pWpM > p(*i+*2). 

(3.8) p(*i)p(*2) = CS^RCS^R > CS tl+t2 R = p^+ t2 K 

But Theorem I3~3l implies that pC^P^) < p(*i+*2) for all large enough ti,i 2 , 
since the concatenation of two C-heavy paths is again a C-heavy path. As are 
periodic, it follows that V [tl) V {t2) < V (tl+t2) for all t 1 ,t 2 , and we obtain the claim 
combining this with the reverse inequality. □ 

Formulas (13.21) and (I3.4j) imply that if A is S- visualized then all rows of R or 
columns of C with indices in the same cyclic class of C coincide. Hence, when 
working with V^we can assume without loss of generality that all cyclic classes 
have just 1 element and consequently, that Sn c n c is a permutation matrix. This 
captures the structure of P«, which form a cyclic group of order 7. 

As usual ei denotes the vector which has all coordinates equal to except for 
the ith which equals 1. For the rows V^) and columns V.^ of we have: 

(3.9) e JvW=7>V, V^e 3 =V { f. 

Next we study the periodicity of Pw i n more detail. It turns out that the 
columns and rows of p( t+1 ) with indices in N c can be obtained from those of 
by means of a permutation on cyclic classes, while the rest of the columns (or 
rows) are max-linear combinations of the critical ones. We start with the following 
observation on the spectral projector V^ ' := CP, which can be found in [TJ[7]. 

Lemma 3.5. P.[ 0) = R. L . and V { ° ] = d for all i £ N c . 

Proof. As (^4 7 )ij = 1 for all i € N c , we obtain efC = Ci. > ef and Pe, = R.i > e,. 
We see that 

p (0) = e fCR> ejR = Ri., 

(3.10) 

p\ 0) = CRe, > Ca = d. 

But CP < ((^4 7 )*) 2 = (^4 7 )*, which implies the reverse inequalities V- 0) < R4. and 
Pi 0) <d. □ 
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Theorem 3.6. Let A e R" x " have A (A) = 1 and be S -visualized. If [i] -> t [j], 

(3 11) •p(*+ s ) _ -p( s ) -p( s ) _ •p(*+«) 

for all s,t > 0. 

Proof. We prove the first equality of (|3.11[) . Using the group law we assume that 
s = 0. We also assume that Sn c n c is a permutation matrix, then efS* — ej. Using 
this and Lemma [3^51 we obtain 



V (0) = e T v (0) = T R = e T 3 t R < 

(3.12) J 1 

<efCS t R = efV^ =Vl 

Analogously we have vfP < "Pj 7 . Multiplying this inequality by and using 

the group law and periodicity, we obtain that < "Pj 7) = vf ] . Combining this 
with (|3.12[) we obtain the desired property. □ 

Corollary 3.7. Let A G R^ xn have X(A) = 1. Then vf = (S'i?)*. and vf = 
(CS*)., for all ieN c . 

Proof. We assume that A is S'-visualized, and we also assume that Sn c n c is a 
permutation matrix. For [i] — > t [j] , Theorem 13.61 and Lemma 13.51 imply that 

V f = efV^ = e Jp<°> = ejR = eJS l R = (S 1 ^., 
Vf = V (t) e, = V { - 0) e l = Ce % = CS'e, = {CS l ). 3 . 
The claim is proved. □ 

Corollary 3.8. Let A G M" x " have A (A) = 1. For each k = 1, . . . ,n there exist 
ctik and Pki, where k G N c , such that 

(3.13) = a ik vf, V® = favW. 

i£N c ieN c 

Proof. By Corollary Mi we have vf = (6*5')., and vf) = (S^i?)*. for all i G N c . 
Eqn. (|3.13p follows directly from = CS t R, the coefficients onu (resp. ft,) 
being taken from the fcth column of R (resp. the A:th row of C). □ 



4. NACHTIGALL EXPANSIONS 

In this section we show that the powers A 4 of A G M™ x ™ can be expanded for 
£ > 3n 2 as sum of CSR products. This establishes a more general algebraic form 
of the Nachtigall expansion studied in [HO [H] ■ 

Let A = {aij) G JR" X ". Define Ai = A (A), and let C x = {N U E X ) be a completely 
reducible subdigraph of C(A). Set Ai := A and K x :— N, where TV = {1, ... , n}. 

The elements of a Nachtigall expansion will be now defined inductively for [i>1. 
Namely, we define K„ := N\ Uf^ 1 N{ and A M = (o{J) G K" x ™ by 

(4.1) af 7 . = ( a «' if ^' G ^' 

I 0, otherwise. 
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Figure 2. Formation of first two terms in a Nachtigall expansion: 
a schematic example. Note that the structure of strongly connected 
components of T>(A) is of no use here. 

Further define A M := X(A^). If A M = then stop, otherwise select a completely 
reducible subdigraph — (N^ l7 E^) of the critical digraph C(A^), and proceed as 
above with fi := fi + 1. 

By the above procedure we define K^, A^, A M and C M = {N^E^) for // = 
1, . . . , m, where m < n is the last number fj, such that A M > 0. 

Denote L := U™ i N n and L = N\L. 

For each \i = l,.,.,m, let 7^ be the cyclicity of C^. Since A((A AI /A AI ) 7f ') = 
A(A M /A M ) = 1, the Kleene star := ((^4 M /A M ) 7 ^)* is finite. Define the matrices 
C7 M = (c£) G Rl xn , i? M - (rg) G R n + Xn and = G M^" by 



(4.2) 





otherwise, 15 | 0, otherwise, 
otherwise 



A schematic example of Nachtigall expansion is given in Figure [2j 

It follows from Proposition ^ . 1 1 that each can be scaled to a — 1 matrix using 

a certain vector denoted here by z M . Note that the sets are pairwise disjoint. 

Defining z G by 



(4.3) 



1, if i G r, 



and letting D := diag(z), we obtain that the matrix A := D~ 1 AD is totally S- 
visualized, meaning that all corresponding matrices are Boolean. 

As can be scaled to be Boolean, the sequences of their max-algebraic powers 
{5* I t > 0}, being powers of Boolean matrices when scaled, are ultimately periodic 
with periods 7^. This periodicity starts at most after the corresponding Wielandt 
numbers (fc M — I) 2 + 1 where is the number of elements in JV M 

We proceed with some notation. Denote fi(i) = fi if i G N^, and fx{i) = +00 
if i G L. Denote by n^ it the set of paths on T>(A) which connect i to j and have 
length t. Denote by n^- t the set of paths P G Hfj.t such that mini £ Ar p = /x, 
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where Np is the set of nodes visited by P. (Note that greater values of fx correspond 
to smaller A M .) The paths in 11^ t will be called fi-heavy, since they are C M -heavy 
(see Sect. [3]) in A^. 

Any path P £ Tlij,t with t > n has at least one cycle. Hence there are no paths 
with length t > n that visit only the nodes in L, for otherwise the subdigraph of 
T>(A) induced by L would contain a cycle and the number of components would be 
more than m. We can express the entries of A t — (ajj) and = {afj ) for t > n 
as follows: 



(4.4) 



Oij = «>(%t) = ©w>(n£. t ), t>n 
m 

f = ©#^), t>n. 



Denote 

(4.5) A/f := C^R,. 

These CSR products are defined from A^ and C M in the same way as 

■pit) 

, see Sect. 

[31 were defined from A and C, and it follows that the sequence is periodic with 
period 7 M . Further denote by the greatest cyclicity of a component in C M and by 
n M the number of nodes in K^. The following is a version of Theorem 1331 for N^. 

Theorem 4.1 (/i-Heavy Paths). Let be such that {St, t > T M } is periodic. 

1. For f > 0, 

(4-6) w(U^ t ) < X^Af^h 

2. For t > T m + 2t m (/i m - 1) , 

(4-7) > A^A/^ 

Proof. We can w.l.o.g. assume that A M = 1, since both (]4.6p and (|4.7|) are stable 
under scalar multiplication of A. After this, the claim follows from Theorem l3.3l □ 



In the theorem above, we can choose T M equal to each other and of the order 
0(n 2 ) for all /i. The main result of this section immediately follows now from 
Theorem [O] and Eqn. (|4.4p. noting that T M + 2r M (n M — 1) < 3n 2 for all //. 

Theorem 4.2 (Nachtigall expansion). Lei A £ K" xn . TTiera /or all t > 3n 2 

m 

(4.8) 4 = A t^ t} - 

In particular, 



m 

(4.9) A* = 0At^*). 

i/=i 
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Figure 3. Formation of first two terms in a Nachtigall expansion 
(upper part) and ultimate expansion (lower part): a schematic 
example. The associated digraph consists of two components, de- 
noted T>i and V 2 - 



5. Ultimate expansion 

In this section we construct a different expansion of A 1 which we call the ultimate 
expansion, in order to describe the ultimate behavior of A*. This expansion is 
related to the Nachtigall expansion of Section @] with the selection rule C M := C(A I1 ). 
The latter expansion will be called the canonical Nachtigall expansion. 

The elements of the ultimate expansion will be labeled by o, since we need to 
distinguish them from those of the canonical Nachtigall expansion. For instance, 
we will write A°^ versus A^, C° versus C M and A° versus A p , etc. 

Let A = (dij) G R" x ™. Dchnc A° = X(A), let C° = (N°,E°) be the critical graph 
of A and denote by M° the set of nodes in all components of T>(A) that contain 
the components of . Set A\ :— A and Kl := N. 

By induction for fx > 2, define K° := N\ Uf^ 1 M° (instead of = N\ Uf N t 
in the case of a Nachtigall expansion), and define A° = (a^°) S R™ x ™ by 



(5.i) < = h> if ^' e ^< 

J 0, otherwise. 
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Define A° := A(A°). If A° = then stop, otherwise let C° = (N°,E°) be the 
critical graph of A^, let M° be the set of all nodes in the components of T>(A) which 
contain the components of C°, and proceed with the above definition for fi := fi+1. 

By the above procedure we define K°, A°, A°, C° = (N°,E°) and M° for 
/i = 1, . . . , m°, where m° < n is the last number /i such that A° ^ 0. Note that 
{A°, /i = 1, . . . ,m°} is the set of m.c.g.m. of all nontrivial components of T>(A), 
and each C° consists of the critical digraphs of (possibly several) such components. 

Let 7° be the cyclicity of C°, and let B°, C°, S° and i?° be defined from 
(A°/A°) 7 * in full analogy with P M , C^, and P M in Sect, g] The matrices S° 
can be again simultaneously scaled to — 1 form. 

Essentially in the new construction we contract by the components of T>(A) in- 
stead of the components of C(A tl ) in the case of the canonical Nachtigall expansion. 
See Figure [3] for a visual comparison. 

Denote by X(i) the m.c.g.m. of the component of T>(A) to which i belongs, and 
let A(i) = if {i} is a trivial component of T>(A). For a path P define A(P) := 
maxigATp X(i) where Np is the set of nodes visited by P. Recall that Ily^ denotes 
the set of paths on T>(A) which connect i to j and have length t. Denote by n^° t 
the set of paths P G Iltj.t such that P contains a node in N° and A(P) = A M . Such 
paths will be called fi-hard. Note that they are C(A°)-heavy with respect to A° . 

Denote 

(5.2) u^-.= c;(s;Yri. 

These CSR products are defined from A° and C° in the same way as were 

defined in Sect. [3] from A and C, and it follows that the sequence {Z4*\ t > 0} has 
the period 7°. Denote by r° the greatest cyclicity of a component in C° and by n° 
the number of nodes in K°. The next result follows from Theorem 13.31 

Theorem 5.1 (^i-Hard paths). Let T° be such that {(S°f, t > T°} is periodic. 

1. Fort> 0, 

(5-3) W (H^) < (A°) 4 (W*)y. 

2. F rt>T; + 2r>°-l) , 

(5-4) «;(D^ t ) > (A°)*(^) y , 

Comparing the constructions above with those of the canonical Nachtigall ex- 
pansion (see Sect. |4] assuming that := C(A^)), we see that C\ is the same as C\ 
and Ai is the same as AJ, however, other components and values may be not the 
same. We next describe relation between them. 

Proposition 5.2. Each v = 1, . . . , m° corresponds to a unique fj, = 1, . . . , m such 
that X° = Ap and all components of C° are also components of . 

Proof. Consider the canonical Nachtigall expansion. Note that A M strictly decrease, 
as at each step of the definition we remove the whole critical digraph. Now pick ar- 
bitrary A°, which is the m.c.g.m. of some component of T>(A). There is a reduction 
step when C M for the first time intersects with a component of T>(A) whose m.c.g.m. 
is A°. Then A M = A°, and C M has to contain all components of C°, precisely as they 
are. This proves the claim. □ 
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Further we renumber A° so that v = fx if A£ = A^ , meaning that the numbering 
of A° is adjusted to that of A M . This defines a subset E of {1, ... , m}, such that 
A^ = A° if and only if // G S. 

Corollary 5.3. 7^ is a multiple 0/7° /or eac/i /ieS. 

Proof. With the new numbering, all components of C° are also components of 
by Proposition 15.21 □ 



Unlike /i-heavy paths, /it-hard paths do not cover the whole path sets Hy^ in 
general. However evidently 

(5.5) (A°X=w(ll!t° t ), ifzeiv; ovjeN^. 

From this and Theorem l5.ll we deduce the following. 
Proposition 5.4. Let A G M" xn . For a// 1 > 3n 2 

(5.6) (A;% = A^«, i/i G A£ or j G 

For the sequel we need to establish some relation between connectivity on T>(A) 
and nonzero entries of Ujp . 

We denote by 7 the l.c.m. of all cyclicities 7° of all components C°. Recall that 
we denote by N P the set of nodes visited by a path P and by A(P) the greatest 
m.c.g.m. of a component visited by P. 

Proposition 5.5. Let i, j ' G N, I > and fj, G S. TTie following are equivalent. 

1. (ujP^^O; 

2. For a/Z i = Z(mod 7 ) smc/i i/iai i > 3n 2 , t/iere is a [i-hard path of length t 
connecting i to j ; 

3. For some t = /(mod 7 ) there exists a \i-hard path of length t connecting i 
to j; 

4. For some t = Z(mod 7 ) there exists a path P of length t such that A(P) = 
A,,.. 



Proof. Implications l.<=>2. and 3<=>1. follow from Theorem 15. II and the periodicity 
of Ltjp . Implications 2.=>3.=>4. are evident. It remains to prove 4.=>3. Let k be 
a node in Np which belongs to M°, and let I be a node in N° (that is, a critical 
node) in the same component of T>(A) as k. There exists a cycle containing both 
k and I. Adjoining 7 copies of this cycle to P we obtain a /i-hard path, whose 
length is congruent to Z(mod 7 ). □ 

Now we establish the ultimate expansion of matrix powers, as an ultimate form 

of the canonical Nachtigall expansion. We will write a(t) = b(t) if a(t) — b(t) for 
all t > t' where t' is an unknown integer, and analogously for inequalities. 

Theorem 5.6 (Ultimate expansion). Let A € R" xn . For all /i£S 
(5-7) K)*= AtZtf>. 

Ln particular, 

1* T 



(5.8) A*£©AtW< t 
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Proof. It suffices to prove (|5.8[) . First note that tijp < A/]u for all fi £ E, since any 

/j,-hard path is a /j,-heavy path. As A — ® A^A^i by the canonical Nachtigall 
expansion, it suffices to prove that 

m rp 

(5-9) ©W<0A^) 

m=i pes 

For all /i, i and j such that (A/)! )ij ^ 0, we will show that either /i £ E and 
(U^)ij = (A/]* )y, or there exists i^eS such that Aj, > A p and (Uv )ij 0. 

Assume that either /i ^ E, or (i e S but (A/^( )y > {U^)ij- Theorem l4.1l implies 
that for all I such that / > 3n 2 and I = t(mod 7^) there exist paths P £ 11^ ; such 

that w(P) — A' We are going to show that these paths are not ^-hard. If 

\x £ E then this is immediate. If /1 £ E then by Corollary [53] 7^ is a multiple of 7°, 

and hence / = t(mod 7°). If P is /^-hard, then w(P) < X^QAjp)^ by Theorem 1 5. 1[ 

which implies (Liji )%j > )»j contradicting our assumptions. Hence P are not ji- 
hard, meaning that for any such path there exists v £ E such that \ v = X(P) > A M . 
Applying Proposition 15.51 we obtain that (Liv)ij with A„ > A M . The claim is 
proved. □ 

If A is irreducible, then the ultimate expansion has only one term, which cor- 
responds to its critical graph C(A). In general, it has several terms (up to n) 
corresponding to the critical graphs of the components of T> (A) (or possibly clus- 
ters of critical graphs of the components with the same m.c.g.m.) Thus, Theorem 
15.61 can be regarded as a generalization of the Cyclicity Theorem, see [17] Theorem 
3.9 or p] Theorem 3.108, which it implies as a special irreducible case. 

6. Computational complexity 

Given A £ R" x ™, we investigate the computational complexity of the following 
problems. 

(PI) For given t, reconstruct all terms A/Jx of a Nachtigall expansion with a pre- 
scribed selection rule for C M . 

(P2) For given t: < t < 7, reconstruct all terms Lijp of the ultimate expansion. 

In (PI) we assume that selecting the subdigraph C M of C(A IJi ) does not take more 
than 0(n 3 ) operations. This holds in particular if C M is an arbitrary cycle of C(A^) 
as in [Tl [21]. 

Problem (PI) is close to the problem considered by Molnarova [12], and Problem 
(P2) is extension of a problem regarded by Sergeev [25]. An 0(n 4 log n) solution 
of these problems is given below. It is based on visualisation, square multiplication 
and permutation of cyclic classes. See Semanci'kova [22] [23] for closely related 
studies in max-min algebra. 

Theorem 6.1. For any A £ R" x ™, problems (PI) and (P2) can be solved in 
0(n A log n) operations. 

Proof. (PI): First we need to compute A^, A^ and for all [i. At each step the 
computation requires no more than 0(n 3 ) operations, based on Karp and Floyd- 
Warshall methods applied to each component of T>(A fl ). The total complexity is no 
more than 0(n 4 ). After this, we find all cyclic classes in each C M , which has total 



1(5 



SERGEI SERGEEV AND HANS SCHNEIDER 



complexity 0(n 2 ), and hence the cyclicities 70- of all components of the graphs C M . 
At this stage we can also find a scaling which leads to a total ^-visualization of A 
(and hence all This relies on Floyd- Warshall method applied to each and 
takes no more than 0(n 3 ) operations in total. 

By Theorem 14.21 A 1 ^ admit Nachtigall expansion for all fi and all t > 3n 2 . The 
rows and columns of A* with indices in are determined at t > 3n 2 only by 
A/"' = C^S^R^, since by construction these rows and columns are zero in all terms 

A/"i* } for v > fj,. This means in particular that these rows and columns become 
periodic after 3n 2 time. By repeated squaring A^, A 2 , At, . . . , we reach a power 
A 7 ^ with r > in 2 , which requires no more than 0(n 3 log n) operations. Now we 
can use Corollary 13.71 identifying C^S^ and S^R^ as submatrices extracted from 
columns, resp. rows, of A^ with indices in N^. By Theorem l3.6l we can obtain S^R^ 
from S^R^ and C M from C^S r ^ by the permutation on cyclic classes determined 
by the remainders r(mod 7 CT ) and (t — r)(mod j a ), for each cyclicity 7<j of a 
component of C M . This takes 0(n 2 ) overrides. Finally we compute A/"* = C M 5* i? M 
(0(n 3 ) operations). We conclude that the total complexity for all /x does not exceed 
0(n 4 log n) operations. 

(P2): It is clear that the computation of all prerequisites for the ultimate ex- 
pansion is done like in the first para of the proof of (PI), and takes no more than 
0(n 4 ) operations. After that we use Proposition 15.41 which means that the critical 
rows and columns in each (A°) r for r > 3n 2 are determined only by IA^\ Hence 

the factors of each 14^ can be computed by matrix squaring of A° followed by 
a permutation on cyclic classes and matrix multiplication, which overall takes no 
more than 0(n 3 logn) operations. We conclude that the total complexity for all \x 
does not exceed 0(n 4 logn) operations. □ 

7. Orbit periodic matrices 

Being motivated by the results of Butkovic et al. [6] on robust matrices, we are 
going to derive necessary and sufficient conditions for orbit periodicity and to show 
that they can be verified in polynomial time. 

A matrix A 6 M™ xn is called orbit periodic if for each y G M™ there exists 
X(y) £ R™ such that A t+ ~<° y — (A(j/)) 7 A l y for all sufficiently large t, where (as 
above)7° is the joint cyclicity (l.c.m.) of the critical graphs of all components of 
27(A). 

A sequence {A t y, t > 0} with the above property will be called ultimately linear 
periodic and X(y) will be called its ultimate growth rate. The same wording will be 
used for the sequences {a\j, t>0. We say that a subsequence {a 1 ^ 1 s , s > 0} has 

ultimate growth rate A, if there exists ay 7^ such that a 1 ^ 1 s = aij\ l+1 s for all 
s starting from a sufficiently large number. 

It may seem more general if in the above definition of linear periodicity we replace 
7 by 7(2/). But using the ultimate expansion (|5.8|) we conclude that {A*y, t > 0} is 

ultimately linear periodic if and only if there exists fx £ E such that A l y = X^U^ y. 

As U^ +1 ■* = ujp for all [i and t, we conclude that the exact period of {A t y, t > 0} 
has to divide 7 . 

The ultimate expansion leads to the following properties of the sequences {a\j , t > 
0}, already known in max algebra O [TBI [20]. As above, A(P) is the largest m.c.g.m. 



CSR EXPANSIONS OF MATRIX POWERS IN MAX ALGEBRA 



17 



of the components of T>(A) visited by P, and Hij,t denotes the set of paths of length 
t connecting i to j. 

Lemma 7.1. For each I: < I < 7 , the subsequence {a^ 1 s , s > 0} is ultimately 
zero or has an ultimate growth rate Ajj(Z). In the latter case, for each such I and 
each sufficiently large t = /(mod 7 ) there exists P £ ^ij.t such that X(P) — Ay (I). 

Proof. The ultimate expansion (|5.8[) implies that for each I there exists jjgS such 
that a 1 ^ 1 s = A^ +7 ° s (W^)ij at sufficiently large s, so the subsequence has growth 
rate Ay(i) := A„. The second part of the statement follows from Proposition 

□ 



Lemma 7.2. If P £ Tlij 1 then {a 1 ^ 1 s , s > 0} has ultimate growth rate at least 
A(P). 

Proof. Using Proposition 15.51 cond. 4, we obtain that (tiff )ij 7^ for /z such that 
A p = A(P), and then a'+ 7 ° s > {X{P)) l+ ^ ' s {U^)ij at sufficiently large s. □ 

Denote by L° the set of nodes in the nontrivial components of T>(A). The next 
statement follows from the Cyclicity Theorem [1] [17] . For the sake of completeness 
we deduce it from the ultimate expansion. 

Lemma 7.3. Ifi,j £ L° belong to the same component ofT>(A), then {a«, t > 0} 
is ultimately linear periodic and its growth rate is the m.c.g.m. of that component. 
Further, for each i £ L° and all t there exist k and I in the same component of 
V(A) such that a\ k ^ and a% ^ 0. 

Proof. All paths which connect i to j are /z-hard, with A p = \{i) = X(j). Using 
Theorem 15. II we obtain that only (U^)ij is nonzero and hence the ultimate expan- 
sion reads af ■ = X^U^)^. For the second part note that i belongs to a cycle with 
nonzero weight. □ 

If i is connected to j by a path, we denote this by i — V j. Observe that if i — > j 
then also k — > I for each k in the same component of T>(A) with i and for each I 
in the same component of T>(A) with j. We also denote i -H- j if both i — >• j and 
j — > i (i.e., if i and j are in the same component of T>(A)). In the next theorem we 
describe, in terms of such relations, when the sequences of columns {A 1 a, t > 0} 
are ultimately linear periodic. 

Proposition 7.4. Let A £ R" xn and j £ L° . The sequence {A t Bj, t > 0} is 
ultimately linear periodic if and only if for all i £ L° , i j implies X(i) < X(j). 

Proof. The "only if" part: Lemma [7.31 implies that for each t there exists k such 
that ajy 7^ 0, and the sequence {a^ S7 , s > 0} has ultimate growth rate X(j). If 
the condition does not hold, there exists a path P leading from i to j such that 
X(P) = X(i) > A(j), and by Lemma [72] there is a subsequence of {aL, I > 0} with 
ultimate growth rate at least X(i). 

The "if" part: If the sequence {A t ej, t > 0} is not ultimately linear periodic, 



then some of its entries by Lemma [7.31 have ultimate growth rate X(j) and there 
is a subsequence of {aL, t > 0}, for some k £ N, which has a different ultimate 
growth rate. Lemma 17. 1 1 implies that this growth rate has to be greater than X(j), 
and must be the m.c.g.m. of a component which has access to j. □ 
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We denote i =4 j and say that i strongly accesses j, if i can be connected to j by 
a path of any length starting from a certain number Tij. For example, i strongly 
accesses j if i — > j and the cyclicities of the components of T>(A) containing i and j 
are coprime. Observe that if i =4 j then also k I for each fcoi and I j. Now 
we show that strong access relations are essential for the ultimate linear periodicity 
of all sequences {^4*(e, © ej), t > 0}. 

Proposition 7.5. Lei A G R" x " and i,j G L°. Suppose that {A t e i , t > 0} and 
{A t ej, t > 0} are ultimately linear periodic. Then {A t (ek © e;), i > 0} is a?so 
ultimately linear periodic for all k O i and Z O j , i/ and onZj/ if either i j , or 
j =4 i, or 6oi/i are false but A(i) = X(j). 

Proof. If X(i) = X(j) then both {A t e k , t > 0} and {A l e u t > 0} for all k o i 
and Z ■<-» j have this growth rate and {A t (ek © e;), t > 0} is ultimately linear 
periodic with this growth rate. So it remains to consider the case X(i) < X(j). 
In this case Proposition 17.41 implies j -/} i, and therefore we have to show that 
{A t (ek © e/), t > 0} are ultimately linear periodic for all k ■<-> i and I ■<-)• j if and 
only if i =4 j- 

The "if" part: For each t: < t < 7 , if there exists m G N and si > such 
that a^ 11 ^ 0, then there exists a path P G II m fc :t+Sl7 o. As /c =4 Z, this path can 
be joined with a path from k to I of length S27 , and we get that o^/^ 1+S2 ' >7 7^ 0. 
Using Lemma [721 we obtain a^ S7 ^ for all sufficiently large s and it dominates 
over since it has larger growth rate. 

The "only if" part: The ultimate linear periodicity of {A t (ek © ej), t > 0}, 
for any k -o- i, implies that supp(A'efc) C supp(^4 t e J ) for all large enough t. By 
Lemma [7.31 there is k <-> i such that (A t )ik 7^ 0, hence also (A 1 )^ ^ 0. As we 
reasoned for any t, it follows that i =4 j. □ 

Theorem 7.6 (Orbit periodicity). A £ E™ x ™ is orbit periodic if and only if the 
following conditions hold for all i,j G L°: 

1. i — > j implies A(i) < A(j), 

2. i/ neither j =4 i nor i =4 j f/iera A(i) = A(j), 

or equivalently if {A t [ei®ej), t > 0} are ultimately linear periodic for alli,j G L°. 

Proof. We need only prove that 1. and 2. are sufficient for orbit periodicity, the 
rest relies on Propositions 17.41 and 17.51 

Let T>(A, y) be the subgraph induced by the set of nodes that have access to 
supp(y) := {i: yi 7^ 0}. Nontrivial components V a of T>(A,y) are ordered by 
relation T> ai ■< T> a2 if i =4 j for some (and hence all) i G T> ai and j G T> a2 . 
Consider the maximal components with respect to this relation, by conditions 1. 
and 2. they must have the same m.c.g.m. and it must be the greatest one. We 
denote this m.c.g.m. by A and show that it is the ultimate growth rate of A l y. 

By Lemma I7TT1 for each I: < I < 7, the subsequence {a^ 81 , s > 0} has a 
certain ultimate growth rate if it is not ultimately zero. We have to show that A is 
the maximal growth rate of {a'^ S7 , s > 0} over k G supp(n), for every fixed ieJV 

and I: < I < 7 . Then it follows that A t+ i°y = X 1 °A t y. 

To avoid trivialities we assume that there exists k G supp(y) and such that 
{ a ifc^ 7 ' 5 > 0} is not ultimately zero. Then for some k G supp(n) there exists 
a path P G Hik.t where t = i(mod 7 ) which visits a nontrivial component T> ai 
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of T>(A, y). If the m.c.g.m. of V ai is A, then by Lemma [7.21 the growth rate of 
{°i^ S7 ' s > 0} is not less than A, hence it must be A and we are done. If the 
m.c.g.m. of D ffl is less than A, then T> ai strongly accesses some component T> (J2 
with m.c.g.m. A, and V a2 accesses a node k' £ supp(y). Due to the strong access 
we can adjust the length of the path from V ai to 2? CT2 if necessary, and we obtain a 
path P' £ Hik',t' where t' = ^(mod 7 ). By Lemma T7. 2 1 we obtain that the growth 
rate of {a'^ 7 , s > 0} is not less than A, hence it must be A. □ 

To assess the computational feasibility of condition 2. in Theorem 17.61 we need 
the following observation which uses the ultimate expansion, see Theorem 15.61 and 
Section [5j 

Theorem 7.7. Let A £ E" x ™ be such that the sequences {A l ei, t > 0} are ulti- 
mately linear periodic for all i £ L° . Then A is orbit periodic if and only if the 
following holds for all /1, v £ S: 

(7.1) A p < A„ |J supple,) C P| supple,). 

Proof. If {A t ei, t > 0} is ultimately linear periodic and i £ L° , then A t ei = 
X^hijpei for A M = X(i), since A M must be the growth rate of A l ei and all other terms 
of the ultimate expansion have different growth rates. Then also A (ei © ej) — 
X t JA^ei © X\,Uv ej where A M = \{i) and A„ = (equivalently, i £ M° and 

j £ M°). If A p < X u then {A* (a © ej), t > 0} is ultimately linear periodic if and 
only if A (ei © ej) — X u tlv £j- This happens if and only if 

(7.2) A M < Aj, => supple,) C supp(wW ej ) V* e M°, j £ M°, 

holds for all fi, v £ S. 

Corollary 13.81 implies that any column of resp. lAu is a max-linear combi- 
nation of columns with indices in ^V°, resp. N°, so that the support of that column 
is the union of supports of certain columns with indices in N°, resp. N°. Hence 
we need to check the support inclusions only for i £ N° and j £ N°. Further, 

Corollary 13.71 implies that the columns of U$ (or resp. U.^) with indices in N° 
(or resp. N°) just permute as t changes, so the inclusions need be verified only for 
t = 1. This shows that (|7.2p is equivalent to 

(7.3) A M < X v => supp(Z/Wei) C supp(WW ej ) £ N°, j £ N°, 

which is equivalent to (17.11) . The claim is proved. □ 

Now we give a polynomial bound on the computational complexity of verifying 
the orbit periodicity of a reducible matrix. 

Theorem 7.8. Let A £ R" x ". Suppose that all components of T>(A) and access 

relations between them are known, and X^ and are computed for all \i £ E. 
Then the orbit periodicity of A can be verified in no more than 0(n 3 ) operations. 

Proof. For all pairs /i, v, we must verify condition 1. of Theorem 17.61 and condition 
(|7.1[) . The first of these conditions is verified for the pairs of components of T>(A) 
and it takes no more than 0(n 2 ), if all access relations between them are known. 
To verify condition (|7.ip . we need to compute all unions UieJV° su PP(^ii e%) and 
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intersections f] ieN o s\ipp(ldjp &i) , which requires 0(n 2 ) operations, and then make 
0(n 2 ) comparisons of Boolean vectors, which requires no more than 0(n 3 ) opera- 
tions. □ 

We combine the results of Theorems 16.11 and 17.81 

Corollary 7.9. Given A 6 R" x ™, it takes no more than 0{n A \ogn) operations to 
verify whether it is orbit periodic or not. 



Mn ax, + 



8. Examples 

All examples in this section will be in the max-plus setting 
{— oo}, © = max, ® = +). 

Example 1. We construct the canonical Nachtigall expansion of A* for 

/-l -7 -6\ 

-1 -5 

-7 -5 -1 

\-6 -4 -3 

A. The maximal cycle mean is Ai 



(M U 



.4 



-4 
-3 

V 



We start with A\ 
has two nodes 1, 2 
We proceed by 
obtaining 

(8.2) 



0, and the component C\ 
thus 



A, = 



V 



-co 

-00 
-00 

-oo 



) and (2,1). 


The cyclicity is 71 := 2. 


in first two 


rows and columns to —00 


—oo — oo 


— oo^ 




— oo — oo 


—00 




— oo —1 


-3 




— oo —3 







The maximum cycle mean is A2 = — 1, and the component Co, has one node 3 and 
one edge (3,3). The cyclicity is 72 = 1. 

Now we set everything to —00 except for the entry (4,4): 

/— 00 —00 —00 — oo\ 

(8.3) A 3 = -°° -°° 

—00 —00 —00 —00 

00 —00 —00 —2 I 

The maximum cycle mean is A3 = —2, and the component C3 has one node 4 and 
one edge (4,4). The cyclicity is 73 = 1. 

We obtain matrices Si, S2 and S3 which correspond, respectively, to Ci, Ci and 

C 3 : 



(8.4) 



Si 



s 3 



—00 


—00 

— 00 

/— 00 

—00 
—00 

— 00 





—00 
—00 
—00 

—00 
—00 
—00 
—00 



-00 
-00 
-00 
-00 

-00 
-00 
-00 
-00 



—00 
—00 
—00 
—00 J 

— 00 \ 

—00 
—00 

J 



So 



—00 


—00 


— 00 


—00^ 


—00 


—00 


—00 


—00 


—00 


—00 





—00 


—00 


—00 


— 00 


—00 J 



These are Boolean matrices in the max-plus setting, with entries 00, instead of 
0, 1. 
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Further we need to compute Kleene stars (A 71 )*, ((A 2 - A 2 ) 72 )* = (1 + A 2 )* 
and ((A 3 - A3) 73 )* = (2 + A3)*, and construct matrices d, Ri, C 2 , R2, C 3 and 
i?3 . The critical parts of the Kleene stars are shown below, the rest of the elements 
being denoted by • as we do not need them: 

(8.5) 



(A 2 )* = 



— 00 \ 

—00 
—00 

/ 



(1® A 2 )* 





-1 - 


-5 -4\ 




-1 


- 


-6 -5 




-5 


-6 






I" 4 


-5 


• • J 




/ • 




—00 


• \ 






—00 




—00 


—00 





-2 


v • 




-2 


• / 



(2® A 3 y 



Further we compute the Nachtigall matrices 
^ (0) = d ® Ri = 



v- 



(8.6) 



-1 
1 

5 -6 
4 -5 

/-l 

-6 
-5 



-0 
-6 
-10 
-9 



-1 

-5 
-4 



y— 00 



-5 



-00 



-00 



-9 

-8/ 
-6 
-5 
-11 
-10 



(- 



M 0) = C 2 ® R 2 



(o) 



C3 ® i?3 



•DC' 
-OO 
-OO 
-OO 

-OO 
-OO 
-OO 
-OO 



-00 
-00 
-00 
-00 

-00 
-00 
-00 
-00 



— OO 
— OO 



-2 

—00 

— OO 

—00 
—00 



-5\ 
-4 
-10 
-9/ 

oo\ 

OO 1 



The Nachtigall expansion starts to work already at t 
-1 -5 -4\ 



A = 



(8.7) 



-1 
-5 
4 



V 

/-i 



- 

-6 
-5 



A 3 = 



V 



-1 

-5 
-4 



-4 

-V 

- 5 \ 

-4 
-5 



A/f 5 8 (-2) ® A/" 2 W © (-4) ® W ; 



-2 

-4/ 
— oo\ 

—00 
—00 

/ 
2. Indeed, 

(o) 



r(o) 



Starting from i = 4 the third term can be forgotten: 



(8.8) 



A 4 = 



/o 


-1 


-5 


-4\ 


-1 





-6 


-5 


-5 


-6 


-4 


-6 


I" 4 


-5 


-6 


"8/ 



-AA 1 (0) ®(-4)®AA 2 (0) . 
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The ultimate periodic behavior starts after T{A) = 10: 



(8.9) 



A w = 





-1 


-5 


-4\ 


-1 





-6 


-5 


-5 


-6 


-10 


-9 


V-4 


-5 


-9 


"8/ 



(0) 



Example 2. The following example will illustrate the ultimate expansion: 



f-2 





-3 


-7 


— oo 


— oo 


— oo\ 





-2 


-5 


-7 


— oo 


— oo 


— oo 


-9 


-7 


-9 


-8 


— oo 


— oo 


— oo 


-9 


-6 


-4 


-4 


— oo 


— oo 


— oo 


-8 


-5 


-5 


-4 


-1 


-7 


-5 


-7 


-8 


-5 


-6 


-3 


-6 


-8 


V-6 


-4 


-9 


-3 


-5 


-5 


-5/ 



(8.10) 



We compute the elements of the ultimate expansion. 

Firstly, A\ = A, Ai = 0, and the component C\ = C(A) consists of two nodes 
1, 2 and two edges (1, 2) and (2, 1). 

On the next step we set all entries in the first four rows and columns of A to — oo, 
thus obtaining matrix A%. Its essential submatrix with finite entries is extracted 
from the remaining rows and columns 5 to 7: 

-1 -7 -5> 
-3 -6 
-5 -5 



(8.11) 



2 ess 



5 to 7 x 5 to 7. 



We compute A2 = — 1, and the component C2 = C{A 2 ) consists of the loop (5, 5). 

The components C\ and C2 determine Boolean (i.e., 0, —00) matrices S\ and S2. 
We compute (A 2 )* and (1 <g> A^)*: 



(8.12) 



(A )* = 



(0 


-2 


-5 


-7 


—00 


—00 


— 00 \ 


-2 





-3 


-7 


—00 


—00 


—00 


-7 


-9 





-12 


—00 


—00 


—00 


-6 


-8 


-8 





—00 


—00 


—00 


-5 


-6 


-6 


-5 





-8 


-6 


-8 


-7 


-8 


-7 


-4 





-8 


V-4 


-6 


-7 


-7 


-6 


-10 


J 




(8.13) (1 ® A2 ess )* = I -2 -(i I . ") to 7 x 5 to 7. 

Next wc build matrices C°, Rl and C% , R% whose essential parts are shown below: 



(8.14) 



(8.15) 



W ess — 
ll l ess 



-2 


V 


-6 


-2 


-9 


-8 


-2 


-5 


-7 


-2 


-3 




°2 ess — 


(0 


-2 


no _ 
"2 ess 


(0 


-6 



1 to 2x 1 to 4. 



-4) , 5 x 5 to 7. 



1 to 7x 1 to 2, 
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Using HOP) and (|8J5]) we compute the ultimate terms U { "> = C{ ® flf , = 
CI® SI® R° and = C% ® R 2 . The ultimate expansion starts to work at t = 9, 
meaning 

( 8 - 16 ) 

A 11 = u[ 1] © (-11) A 1Z = 

In this case, the canonical Nachtigall expansion starts to work already at t = 3, 
and after t = 4 only the first two terms are essential. In this expansion, matrix A 2 
results from setting only the first two (instead of four) columns and rows to — oo. 
The second Nachtigall term J\f 2 is equal to C 2 ® R 2 , where C 2 — Cf and R2 ^ R% 
has essential part 

(8.17) R 2ess = (-A -3 -6 -4), 5 x 3 to 7. 

Example 3. We illustrate the total periodicity. Let 

I —00 1 —00 
— 00 —00 1 

—00 
—00 
—00 
— 00 



AO) 



-io)®wi 0) 

(-12 



(8.18) 



.4 



—00 
1 

— 00 
— 00 



—00 
—00 
-2 
-2 



—00 
— 00 
1 

— 00 
— 00 
— 00 



—00 

— 00 
— 00 
—00 
—00 




— 00 \ 

—00 
—00 
—00 


—00 



(8.19) 



B = 



/— 00 


1 


— 00 


—00 


— 00 


— 00 \ 


—00 


—00 


1 


— 00 


—00 


—00 


—00 


—00 


—00 


1 


— 00 


—00 


1 


—00 


— 00 


— 00 


— 00 


— 00 


— OG 


-2 


—00 


—00 


—00 





y— oo 


—00 


-2 


—00 





—00J 



Note that A and B are almost the same, except for the entries (6,2) and (6,3). In 
both cases the ultimate expansion coincides with the canonical Nachtigall expan- 
sion, and we have two critical components C\ — (Ni,Ei) with Ni = {1,2,3,4}, 
Ei = {(1,2), (2, 3), (3, 4), (4,1)} and C 2 = (N 2 ,E 2 ) with N 2 = {5,6} and E 2 = 
{(5, 6), (6, 5)}. The eigenvalues are Ai = 1 and A2 = 0, and the cyclicities are 
71 = 4 and 72 = 2, their l.c.m. is 7 = 4. Notr that = Ni and M 2 ° = N 2 . 

In both cases the component with Ai does not have access to the component 
with A2 which is smaller, hence condition 1. of Theorem 17.61 is true meaning that 
all columns of A 1 and B l are ultimately periodic. Condition 2. of Theorem 17.61 
holds for A but it does not hold for B. In particular, there are only paths of odd 
length connecting node 6 to node 3. Hence A is orbit periodic and B is not. 

Consider also condition (|7.ip . We need the terms of the ultimate expansion for 
<(mod 7 ) = 1. In each case there are two terms, which we denote by and U 2 , 
resp. and , for the case of A, resp. B. We have 

/— 00 —00 — 00 — 00 — 00 \ 

-00 —00 — 00 —00 —00 

-00 —00 —00 —00 —00 

— 00 —00 —00 —00 —00 

—4 —3 —6 —5 —00 —00 

\-A -3 -6 -5 -00 -00/ 



.20) 
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(8.21) 



U 2 B 



1— oo 





— oo 


— oo 


— oo 


— oo\ 


— oo 


— oo 





— oo 


— oo 


— oo 


— oo 


— oo 


— oo 





— oo 


— oo 





— oo 


— oo 


— oo 


— oo 


— oo 


— OO 


-3 


— oo 


-5 


— oo 


— oo 


V-4 


— oo 


-3 


— oo 


— oo 


— oo j 



.22) 



U 



2 ess 



u 



2 ess 



-OO 







-oo 



5 to 6 x 5 to 6. 

{5, 6} and j e Mf 



Observe that supp(W^ei) C supp^^e^) for all i e Af| 
{1, 2, 3, 4}, but this condition does not hold for B. 

To see the difference between A t x and B t x, take x = [0 — oo — oo — oo — oo 0] T . 
The sequence {A*a;} is ultimately periodic starting from t = 4 with period 4 and 
growth rate 1. In particular the last component of {A*^} yields the following 
number sequence for t > 4: 



(8.23) 



{A t x) 6 = {1 1 1 1 5 5 5 5 9 9 9 9...}, t > 4. 



The sequence {B t x} is not ultimately periodic. In particular, the last component 
of {B t x} yields the following number sequence for t > 2: 

(B*x) 



(8.24) yn x )6 

which can be expressed as 



3.25) 



{0 1 4 5 8 9...}, t>2, 



0, if t is even and t > 2, 
t - 3, if t = 4k + 3 and k > 0, 
t-5, if t = 4k + 5 and fe > 0. 
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